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Abstract: We study the Polyakov line in Yang-Mills matrix models, which include 
the IKKT model of JIB string theory. For the gauge group SU (2) we give the exact 
formulae in the form of integral representations which are convenient for finding the 
asymptotic behaviour. For the SU{N) bosonic models we prove upper bounds which 
decay as a power law at large momentum p. We argue that these capture the full 
asymptotic behaviour. We also indicate how to extend the results to some correlation 
functions of Polyakov lines. 
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1. Introduction 

The proposal [1] for a matrix model describing IIB string theory has triggered consid- 
erable interest in Yang-Mills matrix models. The models are defined in section ^ and 
it is the D = 10 dimensional maximally supersymmetric version which is proposed 
as the model of string theory. When the gauge group is SU{2) it has been known 
for some time how to compute the supersymmetric integrals [2-7]. It was believed 
that the bosonic integrals (in which the fermionic degrees of freedom are omitted) 
would diverge because of the fiat directions in which the matrices all commute, but 
the authors of [8] computed them analytically for SU{2) and numerically for some 
other groups, and realised they can converge as long as D is big enough. We es- 
tablished analytically the convergence criteria for the bosonic SU{N) integrals in [9] 
and for the other gauge groups and supersymmetric integrals in [10], confirming that 
the partition function does exist in the bosonic case when D is big enough (except 
for SU{2) and SU{3), D need only be bigger than 2), and that the supersymmetric 
models exist when D = 4, 6 and 10. For further details, see also [11]. 

For SU {N > 2) it is not known how to compute the integrals exactly, but the 
authors of [12] used the supersymmetry to deform the partition function into a co- 
homological theory in which the integrals can be done. Although the relation of 
this model to the original Yang-Mills model remains unproven, the resulting num- 
bers have been confirmed by very careful numerical calculation for small values of 
N [8,13], and support a conjecture [14] based on D brane dynamics. The method 
has also been extended to the other simple groups, and checked numerically for small 
rank [15,16] and see [17]. Unfortunately, it cannot be extended to observables such 
as the Polyakov line which break supersymmetry. 
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Another important line of research has been into overcoming the considerable 
difficulties in obtaining numerical results for larger values of in order to gain 
understanding of the large N limit. Various simulations of the bosonic SU {N) model 
have been reported [18-22], and the Polyakov line and Wilson loop have been studied 
up to = 768 in [23]. For the supersymmetric model with D — A, the Pfaffian is 
positive, and these models have been studied up to A" = 48 [19,24,25], and with 
particular reference to Polyakov lines and Wilson loops [26]. For D — 6 and D — 10, 
the situation is much more difficult since the Pfaffian is complex. Various ways 
of dealing with this have been tried [27-30]. Meanwhile, supersymmetric random 
surface, and geometrical approaches have been investigated in [31,32], as have mean 
field theory and 1/D expansions in [33-38]. 

The purpose of this paper is to obtain rigorous information about the large p 
behaviour of the Polyakov fine for SU{N) with finite A^. The scene is set in section 
2. In section 3 we derive the exact results for the bosonic and supersymmetric SU{2) 
models; this is not new but our derivation gives an integral representation which 
enables the asymptotic behaviour to be calculated very easily. In section 4 we derive 
upper bounds on the large p behaviour for the Polyakov line in the bosonic theory. 
These decay polynomially, and as we shall indicate, there is good reason to believe 
that the upper bounds capture the true behaviour. We also indicate how to extend 
the results to some correlation functions of Polyakov lines. Section 5 contains a brief 
discussion. 



2. Definition of the models 

The Yang- Mills matrix integral partition function, which is obtained by dimensionally 
reducing the Euclidean SSYM action from D down to zero dimensions, is given by 

2d,g^ J HdX.lldTp^e^p (IT.^ + lT^MKpX,,^p]] , (2-1) 

where we adopt the summation convention for repeated indices. The traceless her- 
mitian matrix fields and ipa (respectively bosonic and fermionic) are in the Lie 
algebra Q of the (compact semi-simple) gauge group G and can be written 

a=l a=l 

where {t", a = 1, . . . ,g} are the generators in the fundamental representation. The 
are ordinary gamma matrices for D Euclidean dimensions. The model possesses 
a gauge symmetry 

[/%[/, ^P^^uHJJ, UeG. (2.3) 
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and an SO{D) symmetry inherited from the original D-dimensional Euclidean sym- 
metry of the SSYM. Although the motivation discussed above leads to a study of 
the D = 10, SU{N) supersymmetric integral, it is useful and illuminating to study 
several different versions of the model. Firstly by suppressing the fermions we get the 
bosonic integrals which we will denote by jV' = (ie there are no super-charges) [8] . 
Secondly the supersymmetric integrals can be written ior D — 3, A, 6, and 10, hav- 
ing J\f — 2{D — 2) super-charges. In principle one can integrate out the fermions to 
obtain 

D 

Zd,g = / \[dX^VD,G{X^.)e-K^{-SD{X)) , (2.4) 

where 

1 ^ 

^^W = -4 ETr[^.,^.]' (2.5) 
/i,i/=i 

and the Pfaffian Vd,g is a homogeneous polynomial of degree \Hg. In this paper we 
will be concerned with Polyakov line correlation functions of the form 

L{P) = < Tr exp(iP"X^) > 

- jWdX^T, exp(iP'^X,)7'^,G(^^)exp i\Y.^ ■ (2-6) 

Using the SO(D) symmetry we can always consider P"^ — (p, 0, . . .) and denote the 
corresponding Polyakov line by L{p). 

It is convenient to express the Lie algebra Q using the Cartan-Weyl basis 

{H\E^}, (2.7) 

where i runs from 1 to the rank r and a denotes a root. In this basis 

[H\W]=Q, [H\E°']= a'E'' (2.8) 

and 

[E", E^] = Naf3E''+(^ if a + /3 is a root 

^2\a\'^a-H ii a = -p (2.9) 
= otherwise. 

Here E~°' — {E°'y, and the normalisation is chosen such that 

TrWW^S'\ Tr£;"£;^ = 2|aP^5"+^, TrH'E'^^0. (2.10) 
We can expand any matrix in the Lie Algebra as 

M = M ■ H + ^M"E" (2.11) 
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with M " = (M")* (but for readability's sake we will only use this expansion when 
necessary). We will denote the set of all roots by A, the set of positive roots by A+, 
and the set of simple positive roots by A^; in addition we will denote a simple root 
by s, retaining a for a generic root. The dual weights Us satisfy 

2^ = Sss'. (2.12) 

The above definitions ( p.7| - p.l2|) fix the normalisation of the long root to be ^/2. In 
the case of su{N), all roots are of equal magnitude \/2. 

Since the integrand and measure are gauge invariant, we can always make a 
gauge transformation ( p.3|) to move into the Cartan subalgebra 

Xd = Xd-H (2.13) 

and reduce the integral over Xd to an integral over its Cartan modes [15] 

f[dX% ^ Qg ( n^AU A^(A^), (2.14) 

a=l \i=l J 

where Qg is the volume of G and the Weyl measure, given by 

Ag(Ad)= n (^D-a)\ (2.15) 

is the generalisation from SU{N) of the Vandermonde determinant factors. 

We will consider mainly the integral ( |2.6| ) without fermions so that Af = and 
there is no Pfaffian. In the Cartan representation the expression ( |2.6| ) takes the form 

Hp) = I ( n^^AU A^(Ab) W rfXfcTr exp(^pAB • H) 
\i=i J k=i 

E^-^^'Zin' + lT.T'lXj.xA (2.16) 



exp 



and since the weights in a given representation correspond to the eigenvalues of the 
Cartan generators, we have 

Tr exp(zpA/) ■ H) = exp^ipXp ■ h), (2-17) 

h 

where {h} are the weights. In the case of su{N), every weight in the fundamental 
representation is a Weyl transformation of every other weight, so we can use Weyl 
invariance to write 

L{p) =N I dXDAUXD)dXk exp UpXo • - ^^fr^^ E I^^T " ^^-i 

^ \ a>0 1*^1 k=l 

(2.18) 

since Ui is the first weight in the fundamental representation. 
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3. Exact calculations for su{2) 



In the simplest case of su{2) it is possible to reduce L{p) to a simple integral rep- 
resentation. This is because there are few enough fields to make good use of the 
rotation symmetry.^ Several authors have used exact calculations of the su{2) parti- 
tion function in various contexts [2-7], and integral representations of the eigenvalue 
density have been given in [40] . In principle one only need take the Fourier transform 
of the densities given in [40] to obtain the Polyakov line but it is hard to extract 
the asymptotic behaviour from these representations. The following calculations 
amount to an alternative representation for the eigenvalue densities from which the 
asymptotic behaviour of L{p) can be easily determined. 

The definitions in section ^ fix the normalisation of the generators, but to avoid 
any ambiguity in the definition (^H) we give them explicitly here; 



(3.1) 



1 /O 1\ 1 fO -A 1 fl 

so that H = (T^ and E"^^ = ^(o"i ± 20-2) • 
3.1 Bosonic 

There is only one positive root and ( p.l8 ) simplifies to 

D-l / D-l 



L{p) = 2 fciXoXlH dXk exp{tpXD/V2) exp I -A^ \xl 

^ k=l \ k=l 

Setting 

Xk = MfcCTi + Vk(J2 + k = l...d, 
where d = D — 1, we get 



Sd-1 



dp'^ 



with 



D-l / 1 1 \ 

K{p) = / dXe'P^ H d^kdukdvkexp -- ^ ^,U,,^, - 4^ - ' 

k=l \ i,j=l J 



where 



U.J 
V 

w 



{u^ + v'^)5ij — {uiUj + ViVj) 

2 2/ \2 
U V — [U ■ V) 



2 I 2 
U + V . 



(3.2) 

(3.3) 
(3.4) 

(3.5) 



(3.6) 
(3.7) 
(3.8) 



^Thc authors of 
matrices. 



also managed to extend this to the integrals over the 3x3 real symmetric 
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Integrating out the ci-dimensional vector ^ gives, up to a constant factor which we 
will fix later, 

K{p) = J dXdudv (det U)'^ exp{ipX/V2 - V/A - X'^W/2) (3.9) 



which, noting that the determinant is det U — W V, we rewrite as 

,^,B,,!5KZ^!M±|MM^^(^,S)^ (3.10) 

A—B2 

where 

J = J dudv d{A - W)d{B - V) (3.11) 

= constB^ e - 5^ , (3.12) 

as one can check by using the rotation invariance. Then integrating out A, and 
scaling Ahy p^, we obtain 

K{P) = P^"^ rdBB'^e-^'' HdAA-'^e-^, (3.13) 

r ( 4 - 1) Jo 

where we have fixed the constant factor by requiring L(0) = 1. Differentiating twice 
with respect to p, we obtain terms which vanish like a polynomial from differentiating 
p^~^ and terms which vanish exponentially from differentiating the integral. Then 
we find 

L(p)^-2^-V^I^|^/-^ (3.14) 

as p — > oo for D > 4.^ 
3.2 Supersymmetric 

For D^4:,6 and 10, the Pfaffian is given by [13] 



where 



-Tr [X^, X,y] [Xj,, Xp] [Xp, Xjj] 



D-2 

2 D-2 



= [6(AV + e,M,,-e,-)] ' , (3.15) 



= 6imV - 2 [v'^UiUm - U ■ v{viUm + UiVm) + U^ViVm) ■ (3.16) 

We note that M and U have the same eigenvectors. Specifically, choosing d — 2 
linearly independent vectors a each perpendicular to u and v, we have 

Ua = Wa, Ma = Va (3.17) 

Uu^v^u-u- vv, Mu^O (3.18) 

Ur = v^v -u-vu, Mv^O. (3.19) 
^For su(2), the partition function diverges when D < 4. 
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The eigenvalues of U on the u-v space are g^-i and qd where qd-iQd = V. 

Proceeding as in the bosonic case and working in the eigenvector basis of U and 
M gives 



D-2 
d-2 \ 2 



^2 

k=l 

d-2 



exp i^-^Wj^a - ^g.-id-i - l^^C -\V- \\'W^ . (3.20) 

As before, we can use the rotation symmetry to do the u and v integrals; the Pfaffian 
terms are pulled down by differentiating with respect to W] and then is integrated 
out to give 

(3.21) 

from which the eigenvalue density p(A) can be read off. Integrating out A gives 

m = -No§^ [ dAA-^-Oe-^y^o (-lj)'^^^e-fl, (3.22) 

where the overall constant factor 

^ 23(3-°)(i.-3)(fl-2)v^ 

is fixed by requiring L(0) = 1. Then, by the saddle point method. 

Lip) ~ (-l)-+^2^--^-Af^p^ exp ( --23p3j (3.24) 

as p — >^ oo. Note that the asymptotic behaviour is exponential and that the power 
law terms present in the bosonic case have disappeared. It is tempting to deduce that 
this is a consequence of the supersymmetry. However it is easy to check that replacing 
the ^{D — 2) power in the expression for the Pfaffian ( p. 151) by any positive integer 
has the same effect of suppressing the power law terms but does not in general come 
from a supersymmetric theory^. The quantity appearing in (|3.15| ) is special in su{2) 
because it is proportional to a derivative of the action when expressed in appropriate 
variables; it is the fact that one of the integrals is then of a total derivative that 
suppresses the power law. 



^We thank Bergfinmir Durhuus for pointing this out to us. 
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4. Upper Bounds on Bosonic Large p Behaviour 

By using the convergence techniques of [9, 10] we can obtain an upper bound on 
the large p behaviour of the Polyakov hne for su{N). We begin by illustrating the 
method on su{2) and then on su{3) before giving the general argument. 

4.1 su{2) 

Starting from ( PT^ ) we integrate out A/5, giving 

, D-l 



/^«^«-»(^-l|)-''(-|^-^- 



(4.l) 



where 

D-l 



|2 

-k\ ! 



j,k=l 

SO that a bound on L{p) is given by 

, -D-l 



k=l 
1 

5i^-i = -^$^Tr [X,,X,]% (4.2) 



1^(^)1 < f / n d^^Q-^ (2 + Q exp (-1^ - 



(4.3) 



In the region of integration which leads to power law behaviour, we must have Q 00 
but Sd-i bounded, so we will eventually need to examine the flat directions of the 
{D — l)-dimensional action Sd-i- First though, in any region in which Q^^ > 
with e a small positive parameter, we have exponential behaviour of L{p), so we 
restrict attention to the region Q^^ < p^^^^*^, giving the bound 

\L(P)\ < ^ jlidXkp-'^'' (2 + ip^^) exp (-^ - So^^'j . (4.4) 
In order to consider the fiat directions, we use the radial variable R defined by 

D-l 

Xk = Rxk, ^ Tr XfcXfc = 1. (4.5) 
fc=i 

Then certainly Q < i?^, so we have the bound 

\m\ < p-'^'- (2 + ip^^) ^ / n (-^ - (4-6) 
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where 

1 

Sd-1 = 5^ Tr [xj, Xkf . (4.7) 
j,k=i 

To estimate the integral at large p we split the integration domain into two: 

7^2: Sd-i<R'-\ (4.8) 

where e is a small positive constant. Since we know that the integrals are convergent 
the R integral in TZi may be evaluated by the saddlepoint method and we get a 
contribution to L{p) which decays exponentially at large p. On the other hand 
the contribution from 7^2 cannot be treated in this way and generates power law 
behaviour. Bounds on the angular integrals in ([4.6|) have been computed in [9] and 
for su{2), they give 



\L{p)\ < p''+'^ { 2 + IpA ^ fdR exp (-^ ) (4.9) 



2"^ J 2 J \ 8i?2 

where e' is another arbitrarily small positive parameter, so that 

\L{p)\ < constp^-^+'" (4.10) 

as long as Z) > 6. 

The difficulty in finding a lower bound is illustrated in this example. In (|4.1| ), 
we have a positive term and a negative term, and each gives the same upper bound 
power law behaviour. One can apply the lower bound methods of [10], and find the 
same power law behaviour as lower bounds for each term (up to arbitrarily small 
parameters e again). These terms originate from the Vandermonde factors and, 
since they each have the same power law behaviour, we cannot rule out cancellation 
between them. 

4.2 su{N > 2) 

Completing the square in ( p.lSj) and changing variables from X to u = X gives, 

f du 1 2 

L{p)= / ~ 1 n ^^fe n {(Q''^v + wQ~^^i) -o^ exp(-pV ■ Q"^ ■ cui) 
J det((5)2 ^^-^ ^g^^ 

X exp I -zy . zy + ^ Tr [X„ Xuf\ , (4.11) 



where 



d 

T^ = Y,\Xkt (4.12) 



k=l 
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The integrand in ( 4.11 ) would be positive definite but for tlie Weyl determinant 
factor which we may bound above 



uji ] ■ a 



aeA4 



< const 



const [u^Q'^u + p'^{uj'[Q~'^uji)~\ 



g—r 
2 



2 



where we have used the relation 



aa^ = const 1. 



(4.13) 



(4.14) 



4.2.1 su{3) 

Now we will temporarily restrict the algebra to su{3) and choose the basis 
si = {J |, —\/\), S2 = (0, V2), in which, abreviating Tj = T^-, 



and 



3{T,T2 + T2T3 + nTi 



4T1 + T2+T3 -v^(r2-T3) 

-V3{T2-n) 3{T2 + n) 

2 4Ti + T2 + T3 



9 (T1T2 + T2T3 + T3T1) 
Since the are positive, we can bound any matrix element by 



^^Q < const\x\\y\ujfQ ^ui, 



(4.15) 



(4.16) 



(4.17) 



which applies in particular to the eigenvalues of Q ^. Taking these bounds together 
with ( [4.131 ), we can bound the expression (|4.11|) for L{p) 

J. D-l {9~r)/2 

\Hp)\< du{u;lQ-'u;,Y/'l[dXk{u;lQ-'u;,Y3-^'^/' Yl a„(z/ ■ i.)^^'-^')/^-" 

k=l n=0 



D-l 



{p^ujlQ~^uJi)''exp{-p^ulQ-^uJi) exp I -z/ ■ z/ + ^ Tr [Xj,Xk 

j,k=i 



(4.18) 



where the a„ are the coefficients of the binomial expansion. We can now integrate 
out u giving some new positive coefficients dn, 

\L{P)\< /n^^'^ E 5„p2"(^fQ^V)"+^/'exp(-j9Vg-V-5^-i(X)). 

(4.19) 



fc=i 
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We now split the integration region into two parts. In the region in which ujJQ^^uji > 
p~2+2e -^[^i^ g small but positive, we have exponential behaviour in p. We therefore 
restrict our attention to the region in which 



u'[Q-'ui < (4.20) 



Writing R as the radial variable = X^X^ + X^X ^ (with summation over /x 
and a implied) we have the bound 



oofQ-'iOi > A, (4.21) 



where A is a constant, since the are positive in ( [4.16|) . Inserting ( [4.20|) and ( [4.21j ) 
into ( [4.19| ) gives 



(9-r)/2 D-l . 2X 

\Lip)\ < J2 -anP^'^ip-^^^r^'" / n ^^'^ ( -"^h ^""p • 

n=0 k=l ^ ^ 

(4.22) 

If < -R < 1, the behaviour is exponential in p. For R> 1, a, bound on the angular 
part of the integral has been calculated in [9], and we can insert the result directly 
into ( [4.22| ) to obtain 

i9-r)/2 / r?\ 

\L{p)\< ^ a„p-^+^(2n+g) / t^i^exp (-A^ j (4.23) 

n=0 ^ ^ 

where e' is another arbitrarily small positive number, and 

kc{D) = 2rD - D-Ar - 5D,^5r,2- (4.24) 

Finally then, we obtain 

\L{p)\ < const p-^-^^~^^~3+'\ (4.25) 

where e" is a new arbitrarily small but positive number. 

These arguments can be extended to some correlation functions of Polyakov lines. 
Consider 

L{p,q) = < Tr exp(ipXi)Tr exp(igXi) > 

= N dXD^G{^D)dXk exp(zpAj) ■ h) exp(igAD ■ h') 

h h' 

x«pf-E%i#i:i^r-so-.V (4.26) 

\ «>o 1*^1 k=l / 
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Clearly the terms with h = h' satisfy the same bound as L{p) with p replaced by 
p + q. After completing the square, the cross terms become 



L{p)= / --7^x11^^'^ n + exp(-a ■ ■ a) 

X exp ( -z. . zy + ^ Tr [X,, X^A , (4.27) 



where, for example (and specializing to su{3)), 

a = puJi — qu!2- (4.28) 
Introducing an arbitrary constant < /5 < 1 we find that 

^^Q- V = ^ (PWm + T2 + T3) + Ti((2p - qr - P\a\') 

+Uip - 2qf - P\a\') + U{p + qf - /^la^) . (4.29) 
It follows that, in the region 

(2p - qf - (3\a\^ = {2p - qf - (3{p^ + pq + q^) > 0, 

(2g - pf - /3\a\^ = (2g - pf - /3(/ + pq + q^) > 0, (4.30) 



we have 



a^g- V > (3\a\\pii + ^2) > const^^, (4.31) 



where /ii^2 are the eigenvalues of Q ^ so that 



\x^Q-'y\ < \x\\y\ J, 2 (4.32) 

and 

{9-r)/2 D-l 
n=0 fc=l 

(4.33) 

where the constants depend on /5. Again, unless cr-^Q^^cr < \af^ we get exponen- 
tial behaviour so proceeding as for L{p) we obtain 

L(p, q) < const\a\-''''^''-'^-'^'"' . (AM) 



Keeping (3 small but fixed, we see that in the regions not satisfying ( 4.30 ) a is almost 
proportional to a root. Since the roots are weights for the adjoint representation 
there is screening and it is easy to see that we just get a constant upper bound 
instead of 
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4.2.2 su{N > 3) 

We now need to generalise the formula [4.16| for u'^Q~^u! to SU{N). First we show 



that for Q defined as in ( [4.12[ ) 



detQ = ^ T^r--Ta., (4.35) 

{ai,---,Or}e5{r) 

where S{n) is the set of all n-tuples of positive roots a which are linearly independent. 
To see this, we use the fact that the determinant is a homogeneous function of the 
Ta of degree r. To find the coefficient of T^^ ■ ■ ■ T^,^ for any choice of the {ai ■ ■ ■ a,,}, 
set = l,'i = l,---,r and = otherwise. Then Q becomes 

Qv = \Y. "^"i (4.36) 

A;=l 

and 

det Q = det(a;i, ■ ■ ■ , a^) det(ai, ■ ■ ■ , , (4.37) 

where (ai, ■ ■ ■ , a^) is the matrix whose columns are the ctj. If {ai, ■ ■ ■ , a^} are linearly 
dependent then det (5 = 0. If {ai, ■ ■ ■ , ar} are linearly independent, then, since for 
su{N) any root is a Weyl transformation of any other root and any positive root can 
be written in terms of simple roots as 

a = Sj + Sj+i H h Sk-i + Sk (4.38) 

with 1 < j < A; < r, we can use a series of Weyl transformations and column 
operations to reduce 

(«!, ■■■,ar) ^ (±Si, ■ ■ ■ , ±Sr). (4.39) 

Any minus signs can be dropped since we are only after the square of the determi- 
nant. The product of (si, ■ ■ ■ ,Sr) and its transpose is just the Cartan matrix whose 
determinant is r + 1 so 

detQ = ^, (4.40) 

which completes the proof of ( [4.35| ). 
Now we show that 

, ,Tn-l, , ^ ^s={ai,-,«r-l}g.S(r-l) ( Iv) '^O'l ' ' "^Qr-l /'/l/IlN 
iO^Lj LOi - Z — , I4.4ij 

where the Ug are non-zero integers. This will allow us to use the same form of bound 
( ^TFD as for SU{3). 
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It is helpful to use a basis in which 



r + 1 



/1\ 



\0/ 



(4.42) 



In this basis 



u^Q a;i = -Q ^/detQ, 

r + 1 



(4.43) 



where Q"*" is the matrix of cofactors. Our task then is to compute the cofactor 
Q'^i I, and as this is a homogeneous function of the of degree r — 1, we can set 
T^jjj = ■ ■ ■ — T/3(^_i) = 1 and = otherwise, as before. Then 



g^ii = 2i-'-detg, 



where the matrix Q is given by 



/lO 





n) 



Vo 



with 



Proceeding as before we note that 



det Q — det 



/I ••• \ 





• Ai) ■ ■ ■ A*--!) 



/ 



So, using linear column operations, we can restore the to to get 



det 



/I ••• \ 





: • • • P(r-l) 



vo 



r + 1 



/ 



(4.44) 



(4.45) 



(4.46) 



(4.47) 



det(a;i,/3(i),---,/3(,_i)). (4.48) 



If the are linearly dependent, then det Q — 0. Using the same procedure as 
before, we obtain 



det(a;i, • • • , P{r-i)) = det{wuji, S(i), • • • , s^r-i)), 



(4.49) 



- 14 - 



where w is a Weyl transformation, and the S(j) are r — 1 of the simple roots. We 
note that this determinant can be zero in some representations of the algebra. For 
example, in the adjoint representation in which the weights are roots. However, we 
now show that it is never zero in the fundamental representation. First, since any 
Weyl transformation is the product of Weyl reflections in roots, we have 

r 

wui = uji — 2 ^^pjSj (4.50) 

i=l 

where the Pi are integers. In terms of the simple roots, ooi is given by 

cui = [rsi + (r - 1)S2 + --- + Sr] (4.51) 

so we find 

q 

det{wuji, ■ ■ ■ , S(,._i)) = det{—^S(r) - 2ns^r), ■ ■ ■ , S(r-i)) 

= ±^det(si,---,s,), (4.52) 
r + 1 

where q is an integer between 1 and r, n is one of the Pi, also an integer, and is a 
non-zero integer. Then we have 

det Q = ^ (4.53) 



r 

and, together with the expression ( [4.35|) for det Q, this completes the result ([4.41|) . 

The upper bound for SU{N) now follows in a manner identical to the proof of 
( [4.41|) . The cofactors of Q take the general form 



Qij= E CsT^.---Ta.-r, (4.54) 

s={oi---a,._i)e>S(r-l) 



where the Cs are positive constants, so that we have again the result ( 4.17 ) 



x^Q < const \x\\y\uj'^Q ^lu. (4.55) 

The remainder of the argument follows without modification from the su{3) case, 
and we have again 

\L{p) \ < constp-''^^^-^^-'+\ (4.56) 
where the function kc is given in |4.24| , and e is arbitrarily small. 



5. Conclusions 

For the bosonic su{2) model the exact calculation shows that asymptotically L{p) ~ 
p"^^^ while the upper bound ( p.l4|) demonstrates that this behaviour arises from the 



fiat directions in the action. For larger gauge groups there is unfortunately no exact 
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calculation but again the flat directions lead to a power law bound on L{p) at large 
p (|4.56|) . In the process of proving the bound we throw away the alternating signs 
coming from the Vandermonde in ( |4.11| ); it is possible that in an exact calculation 
these would lead to the cancellation of the leading power of p. However, as we 
remarked previously, this does not happen in su{2) and so there is no reason to 
suppose that it happens for larger groups. We conclude that L{p) has the asymptotic 
behaviour 

L{p) ~ constp-^^^^-^^-3 (5.1) 

at large p. Our proof of this bound uses certain properties that are special to the 
su{N) Lie algebra but, just as the convergence proofs go through for all gauge groups, 
we expect that very similar results for L{p) will hold for all the compact Lie groups. 
This power law property is different from the exponential asymptotic behaviour given 
by the mean field approximation [34] . The reason is that the mean- field saddle point 
computation misses the dominant region of integration at large p which is essentially 
the end-point of the integration domain rather than the saddle point. 

Although the asymptotic behaviour of the supersymmetric su{2) model is quali- 
tatively different, being exponential rather than power law, it is not clear whether this 
extends to supersymmetric su{N). In the su{2) model the power law is suppressed 
whenever there are fermions in the adjoint representation so this is a feature of 5^(2) 
rather than of supersymmetry. The asymptotic behaviour of the supersymmetric 
su{N) models is an open question. 
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